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ABSTRACT 

Using Baker- Campbell-Hausdorff relations, the squeeze and harmonic- 
oscillator time-displacement operators are given in the form 
exp[<57] exp[a(x 2 )] exp[ft(xd)] exp[j(d) 2 }, where a, j3, 7, and 5 are explic- 
itly determined. Applications are discussed. 
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1 Introduction 



In group theory and quantum optics, the use of Baker-Campbell-Hausdorff relations 
I], 0, to write unitary operators in more useful (often normal-ordered) forms, is a 
well-known technique. For instance, the displacement operator, 

D(a) = exp[<W — a* a] , a = a± + ia 2 = (xo + ipo)/V2 , (1) 

and the squeeze operator, 

S(z) = exp 

can be written as 

D(a) = exp 

and 

S(z) = exp 



-za^a) z "aa 



z = re l< ^ — z\ + iz2 , 



- 2 M 



exp[cW] exp[— a*a] 



^e 10 (tanhr )a t a t 



/ l \(H aTa ) 
V cosh r 



exp 



— —e 1 (tanh r)aa 



exp 



^e 10 (tanh r)a t a t 



(coshr)" 1 / 2 f (seChr , 1)n ( a t)"(o i' 
ZZk n! 



n=0 



x exp 



1 



e 1 (tanh r)aa 



(2) 



(3) 



(4) 



(5) 



However, such transformations are not as well-known in x — p space, in terms of 
the variables 

x = -j=[a + a i ) , d = ip = -^=(a-a)) , (6) 
[a, a f ] = 1 , [x, d] = -1 . (7) 
Granted, the displacement operator is easily expressed as 



D(a) = exp[— ixopo/2] exp[ip x] exp[— x^d] 



(8) 



but the squeeze operator 



S(z) = exp[- Zl (xd + 1/2) + iz 2 (x 2 + d 2 )/2] 



(9) 
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is not. If it were, it could be easily applied to any wave function, using the operator 
properties 

exp[cd]h(x) = h(x + c) (10) 
exp[r(xd)]h(x) = h(xe T ) (11) 

eMc(d 2 )}h(x) = p^£/ x P - ijL ^r- h (y)dy. (12) 

In Section 2 we describe the BCH method to obtain such transformations. In 
Section 3 we explicitly apply it to the squeeze operator, and in Section 4 we do the 
same for the unitary, harmonic-oscillator, time-displacement operator 

T = expH^a + 1/2)] = exp[-i(x 2 + <9 2 )/2] . (13) 

In Section 5 we give examples and indicate further work. 

2 The Method 

There are a number of early papers that deal with this subject, but a complete 
description, that we follow, was given by Wei and Norman [[|. It has been used 
widely. (See, e.g., Refs. [|], [J.) Consider a unitary operator U(t) in terms of the 
exponentiations of /, x 2 , (xd), and (d) 2 times a parameter, t: 

U(t) = exp[t{hl + hx 2 + h{xd) + b 4 (d 2 )}} , (14) 

where the 6« are not functions of x or d. Set this equal to the following ordered 
product: 

Ui(t) = exp[5I] exp[iax 2 } exp[/3(xd)] exp[i7(<9) 2 ] , (15) 

where a(t) (not to be confused with the coherent state label a), (3(t), j(t), and S(t) 
are the functions to be determined. (Note that even in the case where U contains 
only J, x 2 , and (d) 2 terms, one still needs a (xd) term on the right since this operator 
is needed to close the algebra.) 



Take the time derivative of both sides of the equality U = U\, and then multiply 
on the right by W and Jj\, respectively. This yields 

bj + b 2 x 2 + h{xd) + b 4 {d 2 ) 

= [e 25 e _/3 ] {tax 2 + $ exp [iax 2 ] (xd) exp [— iax 2 ] 

+ A f exp[iax 2 } exp[P(xd)](d 2 ) exp[— (3(xd)] exp[— iax 2 } + 5} , (16) 

where "dot" signifies 4. The factor [e 2<5 e _/3 ] is unity by UW = 1. Now the main line 
operators on the right-hand side of the equation should be commuted to the right. 
First do this with 

X = exp[iax 2 ](xd) exp[— iax 2 } . (17) 

Using 



and 



one has 



e A Be~ A = B + [A,B} + [A, [A, B}}/2 + . . . , (1{ 



[ x \ x d} = -2x 2 , [x 2 , [x 2 , xd}} = , (19) 



X = (xd) - i2ax 2 . (20) 



Similarly one has 



Y = exp[(3(xd)}(d 2 )exp[-(3(xd)} 

= d 2 f] ( ~ 2/ ^ )2n = d 2 e- 2 ? . (21) 

n+o n - 



This leaves to calculate 



Z = exp[iax 2 }d 2 exp[— iax 2 } 

= -i2a - iAa(xd) - Aa 2 (x 2 ) + d 2 , (22) 



where the last line comes from direct differentiation. 



Putting Eqs. (EOT), ([210 , and ( 22|) into Eq. ( TO) , one has an equation with coeffi- 



cients multiplying /, x 2 , (xd), and d 2 . Because these are independent variables, that 



means we can write the coefficients multiplying each of these variables as a separate 
equation: 

61 = 7e" 2/3 (2a) + 5 , (23) 
b 2 = ice -i2a(3 + i^e- 2f5 (-Aa 2 ) , (24) 
h = /3 + 7e" 2/3 (4a) , (25) 
6 4 = i 7 e" 2/3 . (26) 

These are four first-order differential equations in four unknowns. They should 
be solved subject to the boundary conditions that the &i(0) = 0. Then set t — 1 and 
one has the unitary operator in product form. (For the time-displacement case, one 
lets t simply remain t, the time.) 

3 The Squeeze Operator 

For the squeeze operator one has, from Eq. @, 

h = -y , b 2 = , b 3 = -zi , 64 = . (27) 



Put this into Eqs. (|23|)-(|2"6|). With linear combinations of Eqs. ( p4|) to one finds 

a = -2z 2 a 2 - 2z x a + ^ 2 /2 . (28) 

The solution that goes to zero at t = is 

. . z 2 sinhrt , A 

where 

<S(i) = cosh rt -\ — -sinhrt . (30) 
Putting this result into a linear combination of Eqs. (B^) and Q25p, 

/3 = — ^1 — 2^ 2 a , (31) 
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or 



Now Eq. (E6|) gives us 



or 



P = -\nS(t) . 



Z 2 2/3 

7 = T ' 



z 2 sinh rt 
7 = = ol . 



Finally, Eq. (p3|) gives us 



2r 



or 



5= --lnS(i) = - . 

2 w 2 

Therefore, setting t = 1, we obtain the squeeze operator, 

iz2 sinh r 



S = s- 1 ' 2 



exp 



iz 2 sinhr. 2 



exp[— (\nS)(xd)} exp 



2r 5 



(9 2 



where 



o = cosh r H sinh r = e cos — he sin — 

r 2 2 



4 The Time-Displacement Operator 

The harmonic-oscillator time-evolution operator can now be calculated in 
way, but with 

h = o, b 2 = -~ , 6 3 = o , 64 = +^ • 

Then the solutions are found in the same way, this yielding, respectively, 

tant 



a = -2a 2 - 1/2 



a(t) 



ap = l/2 + a, /3(t) = -ln(cost) , 
7 = e 2/3 /2, 7 (t) = [tant]/2, 
5= -a, 5(t) = -pn(cost)]/2 . 



This means that the harmonic-oscillator time-displacement operator is 



T = [cos t] x l 2 exp 



— ta.nt(x z 
2 v 



exp [—(In cost)(xd)] exp 



tant(<9 2 



(44) 



This result can be viewed as complementary to others |J in the study of time- 
evolution. 

5 Discussion 



It is enlightening to look at specific examples. 

Using Eqs. (H) and (^) on the harmonic-oscillator ground state, 

^ = vr- 1 / 4 exp[-x 2 /2] , 

one finds 



(45) 



ijj ss = D(a)S(z)ip 

1 exp[-ix p /2] 



where 



exp 



K 



-(x - X f 



Z2 sinh r 
2rs 



1 



2S 2 (1 + i2n) 



ik + ip x 



(46) 



(47) 



This is the most general squeezed state. Setting z to be real and positive yields the 
most commonly studied example: 

(x - x ) 2 



^ = [vr 1/2 r 1/2 exp 



2s 2 



tpox 



where 



s = e . 



(48) 



(49) 



As a test, the time-evolution operator can be applied to the coherent states, which 
are Eq. (|48]) with s — 1. Then one finds 



e -it/2 



~{x — (xq cost + po sint)} 2 



exp[ix(po cos t — xq sin t)] 



exp 



— - (x cos t + po sin t) (p cos t — x sin t) 



(50) 
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Of course, a simpler calculation is possible by replacing a with ae~ lt in the series 
defining the coherent states as an infinite sum of number states. 

But the time-evolution operator can be applied to more complicated systems, for 
instance, the even and odd states of the harmonic oscillator. There, starting from the 
z real and p = squeezed states one can calculate M 



ip s ±(x,t) = Ttp s ±(x) 



(51) 



One finds a closed form expression for ip s ±: 



s cos t — i sin t 



2n l / 2 (l±e- x o cos l ) s 4 cos 2 t + sin 2 £ 



1/2 



exp 



(x — Xq cost)' 



i tant 



s 4 cos 2 1 + sin t 



± exp 



(x + Xq cost) 2 / s 2 — ztant 



— -(tant)x 2 



— j -(tant)x 2 



.(52) 



2 Vs 4 cos 2 t + sin 2 t/ 2 

(Observe that the terms exp[— i(tant)x 2 /2] are necessary to cancel the singularities of 
the terms exp[ix 2 tant/(2 sin 2 1)] when t is an odd multiple of 7r/2.) This then yields 
an analytic description of the probability densities as a function of x, t: 

exp [— (x 2 + xl cos 2 1) j d 2 } 



p s ±(x,t) 



n 1 / 2 d[l±dexp[-x 2 ) /sV 



f 2xxo(cost) \ /2xxosin/f: 
cosh ± cos 



where 



d? 



j2 2 2 j. , • 2 j / 2 

a = s cos t + sin t/s 



V d?s 2 



(53) 



(54) 



Finally, it is amusing to apply the simple time-evolution operator for a particle in 
a box, 



T = exp[+*(<9)72] , 
to a number eigenstate of a particle in a box. One finds 



(55) 



T sin(7rnx) = sin(7rna;) exp[— in n /2] 



(56) 



the correct time-evolution of a number state. Here the operator is repeating the 
calculation for a continuous set of boxes along the real axis. 

These techniques can be applied elsewhere, such as in obtaining the time-evolution 
operator for a system with different potentials. For example, the time-dependent 
system g, | 

V(x, t) = g {2) (t)x 2 + g w (t)x + g {0) (t) (57) 

can be studied. 
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